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from improved staggered quarks * 
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We compare calculations of Bk with improved staggered quarks (HYP, Asqtad) and demonstrate the improved 
scaling behaviour that this gives rise to over previous calculations with unimproved staggered quarks. This 
enables us to perform the calculation of Bk on the MILC dynamical configurations {uf — 2 + 1), for which we 
give preliminary results. 



1. Introduction 

Indirect CP violation in the system of neu- 
tral kaons is parametrized by Sk- Experimen- 
tally it is a very well measured quantity. On 
the other hand, its theoretical expression depends 
on certain CKM matrix elements about which 
we would like to obtain information and on the 
hadronic matrix element between A'° and K''^ 
of the effective hamiltonian describing processes 
with AS = 2, 



rrAS=2 ^ 



with 
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The matrix element of this operator is usually 
normalized by its VIA value, defining Bk as the 
ratio 



\{KO\sj^j5d\0){0\s-f^j5d\K^)' 



(3) 



The theoretical error associated with the cal- 
culation of Bji is the main source of uncertainty 
when one tries to use the experimental value of Ek 
to constraint the CKM matrix, so improvement 
in the determination of this parameter is crucial 
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to gain information about the unitarity triangle. 
The lattice calculations have an important role to 
play in such reduction of errors. 

2. Calculation of using staggered 

fermions 

The staggered quark formulation has the ad- 
vantage of conserving some chiral symmetry at 
nonzero lattice spacing. This residual symmetry 
forbids the mixing with operators of different chi- 
rality leading to a much simpler renormalization 
process. In addition, dynamical calculations with 
staggered quarks are feasible with present com- 
puters. 

The value of Bk used by phenomenologists 
at present in CKM studies was obtained by 
the JLQCD collaboration ^ using unimproved 
staggered quarks in the quenched approxima- 
tion. The value given in P is B^^^{2 GeV) = 
0.628(42). In spite of being considered the bench- 
mark of the calculations of Bk, the determination 
by the JLQCD collaboration has several draw- 
backs. First, it presents large scaling correc- 
tions, so the extrapolation to the continuum is 
not as reliable as it could be with a better scaling 
behaviour. The main uncertainty is quenching, 
whose possible impact was estimated in [2] using 
chiral perturbation theory to be of the order of 
15%. The SU{3) breaking effects are not incor- 
porated in the calculation either, but this prob- 
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Figure 1. Scaling of B^^^{2GeV) with a for im- Figure 2. B'j^^^ as a function of the timeshce 

proved staggered actions compared to the JLQCD (wsea = 0.01/0.05). 

unimproved staggered resuhs. The stars are val- 
ues calculated by renormalizing the unimproved 
bare matrix with the HYP Wilson coefficients. 
See text for a discussion of renormalization. 



ably has a minor effect We show here that 
the large scaling violation can be corrected by 
using improved staggered fermions. On the other 
hand, the dominant error, associated to quench- 
ing, can be reduced by performing dynamical sim- 
ulations. And, finally, the use of kaons made of 
two non-degenerate quarks can account for the 
SU{3) breaking effects. This three-fold improve- 
ment of the staggered determination is the goal 
of our calculation. 

3. Improved staggered actions 

We first check the impact of using improved ac- 
tions in the quenched approximation. We study 
two such actions, the HYP action W and the Asq- 
tad action 4 . Fig. ^ shows our results for unim- 
proved staggered, Asqtad and HYP actions along 
with the JLQCD results, for several values of the 
lattice spacing. 

Two comments are in order in relation with 
Fig. First of all, we match kaon masses at 
a given f3 to those of the JLQCD collaboration 
to make a clear comparison. Secondly, the Wil- 



son coefficients needed in the calculation of the 
renormalized Bk with the Asqtad action are not 
available at present 5 . Hence, as a first estimate, 
we use the Wilson coefficients for HYP 6 in the 
Asqtad case in Fig. ^ This is incorrect, but large 
differences are unlikely because the 1-gluon ver- 
tex is the same for HYP and the Fat7 smearing 
which goes into Asqtad. 

Figfflshows that discretization errors are much 
smaller with the improved actions. This improve- 
ment is largely due to the change in renormaliza- 
tion factors, not to a change in the bare matrix 
elements. The contribution of several matrix ele- 
ments confuses the picture but one way to show 
this is given in Fig. ^ There, the stars correspond 
to the values of S]^^^(2GeV) obtained with the 
unimproved bare matrix elements and the HYP 
renormalization coefficients. 

We expect improved scaling to survive un- 
quenching and therefore allow us to perform reli- 
able dynamical calculations with only a few val- 
ues of the lattice spacing and even obtain useful 
information from a single point simulation. 
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4. Dynamical effects in the calculation of 

Bk 

The next step is to incorporate dynamical ef- 
fects. We performed a dynamical calculation of 
Bk within the Asqtad action, using the configu- 
rations from the MILC collaboration with n/ = 
2 + 1 dynamical flavours 7 . Two ensembles at 
a — 0.125 fm are used: nisea = 0.01/0.05 and 
n^sea ~ 0.02/0.05. From the hadron masses de- 
termined after the simulation, the strange quark 
mass, anis, was found to be 0.04. We make kaons 
on the configurations from degenerate quarks of 
bare mass ms/2 = 0.02. 

Fig. 121 shows an excellent plateau when B'^^'^ 
is plotted as a function of the lattice timeslice. 

In the absence of the correct renormaliza- 
tion coefficients, we again estimate B^^^ {2GeV) 
from the HYP Wilson coefficients for unimproved 
glue. The results we obtained in this way for the 
two different values of the light quark mass, in- 
cluding only statistical errors, together with the 
result quoted by the JLQCD collaboration and its 
corresponding uncertainty, are depicted in Fig. 21 

From the results in Fig. Oit seems likely that 
Bk from full QCD will be close to the quenched 
result. There is also some sign that B^^^ {2GeV) 
is falling with the dynamical quark mass. 

5. Summary 

We found that the large scaling corrections that 
affected previous staggered calculations of Bk are 
reduced by using improved actions (HYP, Asq- 
tad) in the quenched approximation. This im- 
proved behaviour allows us to perform a reliable 
dynamical calculation of this non-perturbative 
parameter on the MILC configurations. We give 
a preliminary estimate for B^^^{2GcV) in Fig. 

01 

In order to get a final dynamical value for 
B^^^{2GeV) and to clear up the role of dynam- 
ical effects in this quantity, we must redo our cal- 
culation using the correct renormalization coef- 
ficients, obtain results for another lattice spac- 
ing and eventually, although its effect is less im- 
portant, consider physical kaons instead of kaons 
made up of two degenerate quarks with ms/2. 
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Figure 3. An estimate of Bj^^^(2GeV) from dy- 
namical lattice QCD. The x axis corresponds to 
the light sea quark mass divided by the (real) 
strange quark mass. The lines represent the 
quenched result from Pp. 
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